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Abstract 

A new exponentially convergent algorithm is proposed for an abstract the first order differential 
equation with unbounded operator coefficient possessing a variable domain. The algorithm is based on a 
generalization of the Duhamel integral for vector-valued functions. This technique translates the initial 
problem to a system of integral equations. Then the system is approximated with exponential accuracy. 
The theoretical results are illustrated by examples associated with the heat transfer boundary value 
problems. 
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1 Introduction 

This paper is devoted to special class of differential equations, which are associated with the first order 
differential equation in a Banach space X 

^-+A(t)u(t) = f(t),u(Q) = u . (1) 

Here t is a real variable, the unknown function u(t) and the given function f(t) take values in X, and A(t) 
is a given function whose values are densely defined, closed linear operators in X with domains D(A,t) 
depending on the parameter t. Equations of the type ([I]) are called abstract differential equations with an 
unbounded operator coefficient possessing a variable domain and can be considered as metamodels of initial 
boundary value problems for parabolic equations with time-depended boundary conditions. 

The variable domain of an operator in some cases can be described by a separate equation, then we have 
an abstract problem of the kind 

^p- = A(t)u(t), 0<s<t<T, 

L(t)u(t) = <S>(t)u(t) + f{t), < s < t < T, ( 2 ) 
u(s) = u Q 

instead of (|TJ) . Here L(t) and are some linear operators defined on the boundary of the spatial domain 
and the second equation represent an abstract model of the time-dependent boundary condition. An existence 
and uniqueness result for this problem was proved in [2J. 

The literature concerning discretization of such problems in abstract setting is rather not voluminous 
(see e.g. [TJ, where the Euler difference approximation of the first accuracy order for problem (|TJ) with 
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the time-dependent domain was considered, and the references therein). It is clear that the discretization 
(with respect to t) is more complicated then in the case of a t- independent domain D(A) since the inclusion 
Vk = y(tk) £ D(A,tk) of the approximate solution yf. at each discretization point should be additionally 
checked and guaranteed. The using of Duhamel integral was proposed in [llj for the problem 



where u(t) is the unknown function with values in a Banach space X, f(t) is a given measurable function, 
A(t) : D{A) € X is a densely defined, closed linear operator in X with a time-independent domain D(A), 
g(t) is a given function with values in some other Banach space Y and d\, do(t) are linear operators. 



g : (0,T) ->■ Y is from L q (0, T; Y), and d x : D(A) -> Y (independent oft!), <9 (0 : D{A) -> Y (can depend on 
i!). Duhamel-like technique allows to transform the problem (J3]) to a system of integral equations possessing 
operator coefficients with ^-independent domains which can be efficiently approximate. It was constructed 
a discretization of high accuracy order for the case when A(t) is a constant operator (i.e. A(t) = A) for the 
problem ([3]) in [llj using this approach. 

The second equation above is just the time-dependent boundary condition with appropriate operators 
d\ , do acting on the boundary of the spatial domain. For this reason we call this equation an abstract 
(time-dependent) boundary condition. Including the boundary condition into the definition of the operator 
coefficient in the first equation we get a problem of the type <m with a variable domain. 

In the present paper we consider the problem ([3]) and build a new algorithm for approximate solution 
which rejects a limitation on the structure of the operator do(t) used in [llj and can be suitable for non- 
constant operator A(t). 

The paper is organized as follows. In Section[5]we transform problem @ to a system of abstract boundary 
integral equations using the Duhamel-like integral. In the Section [3J we construct a numerical method using 
the Tchebychev interpolation for involved unknown functions and the collocation. The main theorem of this 
section shows an almost (i.e. within to a polynomial factor) exponential convergence of the discretization 
for analytical input data. In Section [4] we represent some computational experiment for our algorithm. 

2 Duhamel-like technique for the first order differential equations 
in Banach space 

For the problem ([3J let us choose a mesh ujk = {ti = I * t, I = l,...,K, r = ■£} of K various points on 
[0, T]. Then one can rewrite the problem (j3J) in the following form: 




(3) 



d lU (t) + d Q (t)u(t) =g(t), 
u{0) = Uq, 



Here u : (0,T) -> D(A) C X, f : (0,T) ->• X from L q (0,T;X) with the norm 





+ A(t)ui(t) = /(t), te(ti_i,tj], 



d x ui{t) + d (t) Ul (t) = g(t), 

Ul(tl-l) = 

1 = 1,2,. ..,K, 



(4) 



and 



u{t) = ui{t), te[ti-x,ti}. 

We transform each interval \ti-\,t{\ to the [—1, 1] by change of variables 



t= r_ s+ t 1 + U- 1= r_ {s + 2l _ i)=Msl 
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then wc obtain 



dui{ ^ {s)) + ^ A (Ms))MMs)) = (Ms)), s e (-1, l] , 

diui(4>i(s)) + do(ipi(s))ui(ipi(s)) = g(ipi(s)), (5) 

U J (^(-l))=«i-l(Vi-l(l)) 1 

1 = 1.2,. ...K. 



Let us introduce the following notations: 



vi(s) = ui(i/)i(s)), fi(s) = f(ipi(s)), 
doj = <9oOz(s)), gi(s) = g(M s ))- 



Further let us chose a mesh on segment [—1, 1] with the Chebyshev-Gauss-Lobatto nodes ujn = {s k 
iE^JltE f k = 0, N}. It is well known that max k {6 k } = jr, k = s& — Sfc_i. The problem (JSJ is equivalent 



to 



dMs) = d ,iVi + gi(s), se [-1,1], (6) 

«,(-!) =«l-l(l), 



where 

M,k = Ai(s k ). 

On each subinterval (s^-i, s&] we define the operator ylj fc with ^-independent domain by 

D(A$) = {u e D(A) : d lU = 0}, 
Ag^ = 4,*«VueD(Ag) 
and the operator Bj^ : F — s> -D(A) by 

A Lk (B ltk y) = 0, 
oiBi^y = y- 

For all s £ [—1, 1] wc define the operators 

AW(S) = A% ) se(8 k -i,8 k ], 

Bi(s) = B l>k , s e (s fc _i,s fe ], Vfc = 1, ...,iV. 
Further, we accept the following hypotheses: 

(Bl) We suppose the operator A^ 2 \s) to be strongly positive, i.e. there exists a positive constant Mr 
independent of s such that on the rays and outside a sector Eg = {z 6 C : < arg(z) < 8,9 € (0, 7r/2)} the 
following resolvent estimate holds 

Kzi-A™W-^<T^. (io) 

This assumption implies that there exists positive constants c, K such that [3], p. 103 

||[A( 2 )(s)] K e- AA<2> ( s )|j < cA- K ,A > 0,k > 0. (11) 
(B2) There exists a real positive u) such that 

|| e -AA< 2 >(s)|| < e -^A VA)S g (12) 
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(see [5], Corollary 3.8, p. 12, for corresponding assumptions on A(s) ). 
We also assume that the following conditions hold: 
(B3) 

\\[A( 2 \t)-A^(s)][A^(t)}-^\\ <c\t-s\ Vt,s, 0<7<1; (13) 

(B4) 

\\[A( 2 \t)f[A^(s)}-^ -I\\ <c\t-s\ Vt.aGt-1,1]. (14) 

(B5) 

||d ]]<c. (15) 

(B6) It holds that 



\\[A^m^e- Al2) ^-^B( V )\\^ x dX 



< cVt, r, e [-1,1], < 7. (16) 



Following |llj one can wright down using the Duhamel's technique 

«,(«) = e- A ^ s - s ^v l { Sk _ l )+ T - f e- A ^ {s - x) {-[MX) -A I , fc ] Uz (A)+/ I (A)}dA 

2 "'sfc-l 
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A ( 2 fc ) e- A ^5 ( ^ A) i? i , fe {-a , ; (A)^(A) + , 9 (A)}dA, 



Sfc-l 



(17) 



+ / S e- A -* (s - A) {-[A i (A)-^ fc ] Vi (A) + / i (A)}dA 

+ M*)£ f A (2 fe ) e-^5 (s - A) i? ; , fc {-9 , ; (A> i (A)+ff(A)}dA, 

2 Js k -i 

s e [s fc _i,sfc], fc = 1, ...,N, 
«i(-1) = «i_i(1). 

It was proved in [TT] that under the assumptions Bl- B6 the system (|17D possesses a unique solution in 

y. 

3 Numerical algorithm 

We use the interpolation on the Chebyshev-Gauss-Lobatto nodes in order to construct a discrete approxi- 
mation of ©, CE7]). Let 

N 



Pn(s; vi) = P N vi = ^2 v i( s ]) l j,n(s), 



PN{s;d 0i ivi) = P N (do,i.vi) = >^d j(sj)vi(sj)L jtN (s), 

3=0 

be the interpolation polynomials for i»;(s), 9o,;(s)w;(s) on the mesh ojjv, x = (xq, x\, xjy), Xi G X, and 
V = (2/0,2/1, —,Vn), Vi^Y given vectors and 

AT 

Pjv(s; y) =P N y = J2 VjLjA*) (19) 
the polynomial that interpolates y, where 



Lj N — — ikm^i J=0 _ N 
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are the Lagrange fundamental polynomials. Substituting Pn(t);x) for Vi(r)), x k for Vi(sk), PNivw) for 
9o,i(v) v i(v)i Uk for 9o,z(sfe)f;(sfe) and then collocating in the points s — s k in (jTTj) we arrive at the following 
system of linear equations with respect to the unknowns Xk , y k '■ 



N 



N 



,(0 



X fe ~~ ^ X fc- 



3=0 3=0 
' A- 

3=0 3=0 

which represents our algorithm. Here we use the notations 



V ( k = d o.i( s k) 



(20) 



Sfc-1 



{Aj,* - Ai(rj)} Lj tN (r))drj, 
Bi, k L j<N (r))dr), 

Afle- A ^ s ^B l , k g l {r } )dr 1 - 



(21) 



and suppose that we have an algorithm to compute these coefficients. 

Remark 3.1 Under the assumption that f(t),g(t) are polynomials the calculation of the operators a k j and 
the elements <j) k can be reduced to the calculation of integrals of the kind I s — f* k e~ A<k >f - tk ~^A s d\, which 



can be found by a simple recurrence algorithm: Ii = —I 



Il-i 



A 



(2) 



(t{l-t{_ ie -^)j 



1,2, ...,s, I Q = 



A 



(2) 



I — e Ak Tk ) , where the operator exponentials can be computed by the exponentially 



convergent algorithm from JtJjj, 



After separating of Xq and y^ in (1201) (that we assume are known from the previous step) we have 



N N 

.(0 _ a -A?l %e k (i) , (0 , a (0 , (0 ,v^o .,(0 , j,(0 

c fc -c l - k2 x k _ 1 +a ko x + p k ay + ^a k3 x } + ^PkjVj + <P k , 

3=1 3=1 



= 9a,i{s k ) 



N 



N 



3=1 



3 = 1 



(22) 



fc — 1, iV; Xq — x^ N , 2/q •* — , 



For errors = (z^\, 4'jv)> 4 = (4?i> — ' 43v0 with 4?fc = u K s fc) and zfy = d Q j(s k )vi(s k )-y k 
we have relations 



0-1) 



4!l = 5 o,z(s fc ) 



k = l,...,N; 



N N 



3 = 1 



^o^'+^o^+e-^^^-i+E 



y ('-i) 



3=1 

AT AT 

(0 



^3+E^4!Wfc 

3=1 3=1 



(23) 
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where 



v k 2 



e -A hhi ( Sk -v) {Ai _ Alijl)} {vi{jj ) - P N ( m v l )}dr, 



Sk-l 
Sk 



A^e- A '^ (sk - v) B ltk [doAv)vi(v) ~ P tf (m8o,iv l )]dn. 



(24) 



In order to represent algorithm (|22p in a block-matrix form we introduce the following matrix and vectors: 



/ E x 





V o 





Ex 







4< 2 > T fl 



\ 




E x J 



\ 



p(0 

2/ 



/ 



with i?x being the identity operator in X, the matrix = {&k,j}kj=i with ak.j = [^flV a k,j i^ij ] _7 > 
the matrix = {Pk,j} kj=i w ^ tn @ k ,3 = [^ikV^kd an d the vectors 



5«) 



(25) 



Besides, we introduce the matrix 





-§! £ x o 

V o o o 

with g i _ 1 - e-^^^'nAf^l-Tir 2, iV. 



\ 





-sn-i Ex ) 



(26) 



G 



It is easy to see that for the (left) inverse 



it holds 



/ E x 












\ 


s i 




Ex 




n 

u 


n 
u 


S2S1 




h 










\sn-i ■■■si 


Sn-i ■ ■ 


s 2 • • 


• sn-i 


Ex) 






(Ex 


• 


^ 











E x ■ 











V 


• 


•• E x J 





(27) 



(28) 



Remark 3.2 Using results of JS[ [^J one can get a parallel and sparse approximations with an exponential 
convergence rate of the operator exponentials contained in (S^) -1 and as a consequence a parallel and sparse 
approximation of S . 

We multiply the first equation in (j22|) and the first equation in (|23l) by [Aj 2]' y to obtain a solution of 
and estimating of error. Then, (|22l) . (|23[) can be written in the matrix form as follows: 



(29) 



where 



5040 = CW40 + D d) z (i) + F (d z (j-v + if + ^w, 



A = diag 



5 0)i (si)[4 2 i 3 ]" 7 5 - 5 ^(sJv)[A^]" 7 



(2)- 



The systems (j29j) and (p0| are equivalent to the following ones: 



(30) 



yW = A [(7 — + C«z« + £)(0y(Q + F^x^ + F^y^ + /« 



(31) 



For a vector v = («i,U2, ■■■,vn) and a block operator matrix A = {ajj}|j- =1 we introduce the vector 



norm 



and the consistent matrix norm 



IIMII = llkllloo = max ||u fc || 

Kk<N 



N 



\A\\ 



Ki<N 



En* 



(32) 



(33) 



3=1 



For further analysis we need the following auxiliary result. 
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Lemma 3.3 Under assumptions Bl- B6 the following estimates hold true 

\\\{SV)- x \\\<cN, 
\\\C^\\\<c(^N^nnN, 

|||£> (0 1| | < c T -~N~ llq In N, 1/p +l/q = l, 
ll|A|||<c. 

with a positive constant c independent of N. 
Proof. Due to (B4) we have 

[]|[^-S]^[^S-x]^llI = [[[[^ < f fe ) ]^[^ C f2- 1 ]-' y - ^ + ^111 < 1 -h c|*9 fe . 
Using this estimate, max/j 9k < and (B2) we get further if J < 1 

iiK^nii < i + e -* *a + c ii) + . . . + [ e -**(i + c—r- 1 

<l + (l + C l) + ... + (l + C l) Ar - 1 <^<c7V. 
Using B2, B3 for C (i) we have 

ll|c w lll< f> fci || 



(34) 



(35) 



AT 

r , 



Kk<N ^— ' 2 
j=l 



s fc-l 



AT 



^ L^E / IPSF^^^II IPS -MV)U^\\ \\[A$V[4$]-i\\ \L jtN (r,)\dv 
< \ max / (~(s fc - ?7)) _7 ^(s fc - rj)c^2\L j>N (r])\dri 

— - J s k-l j — l 

<cAjv(J) 2 ^ max ['" (s k - r,) 1 ^ dn 
y 2' i<k<Nj Ski y 

n 

A„ = max } \Lj. n {t)\ 

is the Lebesgue constant related to the Chebyshev-Gauss-Lobatto interpolation nodes. For the matrix D^> 
we have from ([TB| 

N 

{■ A- AT- 



where 



HI^IH ^^T< x avE UM 

„ Sfc AV 

- \ / \\i A $^ e ~ A ®* l ' k ~" )B S$\\ E i^wi^ 

— — J s k-l " — 1 
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< c^A N 



The last estimate is a simple consequence of assumptions (Bl) and (B5). The lemma is proved. 
From the second equation in pip one can write down 



where 



[I - AI)W]yW = A [/- S {i) + C (l) ] + A*W, 



If exists [7- AD^]- 1 we have 

yW - [7 - A£»«]- X A [7 - S« + 5« + [7 - AT^^A^. 

Otherwise one can choose appropriate r so that ||A7?^|| < 1 and in this case it means that there exists 
operator-matrix [7 — ATH^] -1 . Substituting this expression into the first equation in (|3ip we have 



where 



G (0 f W =g (0$(0 ) 

G (D = g«) _ C d) _ D d) {I _ ad(0]-i A [j _ §m + c % 

q(0 =£)(')[/_ ATJ(')]- 1 A + 7 x . 



Similarly one can obtain from ([30 



z(') = [i-ad^aIi-SW + cW 



4° + [7- AD^^A^, 



where 



Lemma 3.4 Under assumptions of Lemma \3.3\ there exists (G^')^ 1 and it holds 

IIKGW)- 1 !!! <cjv, 
IIIQ w lll<c 

tuif/j some constant independent on N . 

Proof. We represent = S^[I X - G^} and estimate \\\G^\\\ with 

G {i) = {s (i) r i c (D + (5(0)-iD(0[/ y _ A7JW]- 1 A(7 X - S« + C<'>). 
We have in the case when exists [7 — AD^]" 1 (this can be always achieved, see comments above) 

lll^lll < \\\(S®r l \\\ ■ |||C (i) ||| + \\\(SW)- l \\\ ■ IIP (0 IIM||A||| • (\\\I X - S^\\\ + |||C«|||) 
and now Lemma 13.31 implies 



(36) 



G^IH < cln^ 



iV 1 ^ V2 



1 T 



TV 1 /?-! 2 



(37) 



This estimate guarantees the existence of the bounded inverse operator (Ix — G±) 1 by the appropriate 
choose of t (to provide |||G^||| < 1) which together with the estimate 1 1| (S^)^ 1 1| | < cN proves the first 
assertion of the lemma. The second assertion is evident. The proof is complete. ■ 
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This lemma and representations of x^ l \ y^ l \ z x and Zy imply the following stability estimates: 

III^IH <cN\\\$W\\\, 
II \zJ l) II I < cN\\$W\\ |. 

Illy (0 lll<civ|||$«|||, 
III^IH <ciV|||^)|||. 

Let IT at be the set of all polynomials in t with vector coefficients of degree less or equal then N. In 
complete analogy with [TJ |9l [10] the following Lebesgue inequality for vector- valued functions can be proved 

Ktj)-Pjv(»7;u)||c[-i,i]= max Jufa) - P N (r);u)\\ < {l + A N )E N (u) (40) 

with the error of the best approximation of u by polynomials of degree not greater then N 

E N (u)= inf max \\u(rj) - p(r})\\. (41) 

Now, we can estimate the error of our algorithm for I's stage. 



Theorem 3.5 Let the assumptions of Lemma Iff.ffl with 7 < 1 hold, then there exists a positive constant c 
such that 

1. For N, K large enough it holds 

1114° II I < cU^} 2 ^ N^lnNENiiA^pvt) + ^-^'InJV^ffiJ,)} , 
1 1|4° II I < c { (J) 2 7 ^ 1_7 ln ^Jv(Ko] 7 ^) + ^iV 1 - 1 /" In NE N (dvi) 



where vi is the solution of (|17p iirci/i i/ie initial condition 
2. TTie equation for x"> can be written in the form 

5 (0 = G (') £ (0 + [5(0]-i Q (0$(0 (43) 
and can be solved by the fixed point iteration 

2(L+i ) =G ! i jcj2) + [5W]- 1 QW$W, fc = 0,l,...;s[J ) -arW*rory (44) 
im'i/l i/ie convergence rate of an geometrical progression with the denominator 

/or JV, if large enough. 
Proof. For z$ we have to estimate in (|38p. 



1 /t\2-7 1 r , 



!>iP |||= max 



2 



Sk-l 



2 



l,k 

Sfc-l 
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< max 

Kk<N 



2 



max 

Kk<N 



" [A (2 fc ) ] 7+1 e" A,( -^ (Sfc ~" ) Bz, fc [^,,(77)i; i (7 7 ) -Pjv^;^)]^ 



Sfc-l 



- c( 2 } i<l<W {sk ~ v) dv ll[Ao] 7 ^(-)-^(-;[Ao] 7 ^)llc[-i,i] 

where dvi(r)) = doj(r))vi(T)). Further, using (|40l) we obtain 

lll^lll < c{^) 2 ^ N 2 ^ XnNEN^A^m) + c^N-^nnNEMidvi). (45) 

The first assertion of the theorem follows from and the second one from (|39p . ■ 
Let us estimate the full error of approximation in the collocating points. 

|||z w |||= max \\vi(s k ) - x['\\ < max \\v[(s k ) - vi(s k )\\ + max \\vi(s k ) - xY> ||, 

l<fc<JV K l<fe<JV l<fc<W K 

where the first summand in the right-hand side of inequality is the error cosed by approximation of the 
initial data for the I's step of our algorithm. Let z$p = [vi{s\] — U;(si)), . . . , w;(sjv) — vi(sn)- Therefore 

< Ill4°lll + lll^lll < e^lll^UI + Ill4°lll = e^HI^HI + |||4°lll 
<---<Ee ( ^>1ll4^l!l<Ee ( ^>Hll4 j) lll 

3=1 3=1 

The same is valid for the error |||4 III = maxi< fc <Ar \\d j((s k )vi(s k ) - y%'\\ 

iii4°iii<Ee ( ^ )CT iii4 J) iii- 

3=1 

Let us introduce the following notation z = (z^\ . . . , z^), zg — (z g , . . . , z g , ). Now we can formulate the 
main result. 



Theorem 3.6 Let the assumptions of theorem \3.5\ hold, then there exists a positive constant c such that for 
N, K large enough it holds 

x " ' iu.iv H/N{is±l,0l ■ vi) t -iv ' ' nut rj N \uvi) f , 

(46) 



4 lll <c|(^) 7 N^lnNENdA^rvt) + ^N 1 - 1 ^\nNE N (dv l ) 
■ | (~) ^ 7 N 1 ^ 1 ' In AT£jv([Ao] 7 wO + -iV 1 - 1 /* In NE N {dv t ) 



1114° lli^ 

where vi is the solution of (jTTJ) wra£/i i/ie initial condition «;_i(l). 
Proof. From the estimates for HIz^HI and |||4 III we have 



K K f 2-7 

< E ei " CT Hl4 J) lll < e c E c (£) ^^IniV^d^-o] 7 ^) + liV 1 - 1 /^^^) 

3=1 3=1 L 



K ( 1- "I 

- C E5i(D 7 ^ 71n ^([^,o] 7 ^) + ^ 1/9 ln^(^) 

3=1 ^ J 

The same is valid for zg 

3 = 1 L J 



(47) 
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4 Numerical example 



In this section we show that the algorithm (|2"0)) possesses the exponential convergence with respect to the 
temporal discretization parameter n predicted by Theorem 13.61 In order to eliminate the influence of other 
errors (the spatial error, the error of approximation of the operator exponential and of the integrals in (1211) ) 
we calculate the coefficients of the algorithm (|20p exactly using the computer algebra tool Maple. 
A special example of the problem from the class (U} is 

du d 2 u 

at = d^ + j{x,t) ' 

u( ,t) = 0, d -^ + b(t)u { l,t) = 9 (t), ^ 
u(x,0) = u (x), 

where the operator A : D(A) £l-)I. X = L q (0, 1) is defined by 

D(A) = {v€W*(0,l):v(0) = 0}, 

Q2 V (50) 

ox 1 

the operators d\ : D(A) — > Y, and d {t) : D{A) — > Y, Y = IR are defined by 

du(x,t) 

dlU = -dx~ ^ ' (51) 
d (t)u = b(t) ■ u(x,t)\ x =i 

and g(t) e L q {0,T;Y) = L q (0,T). 

As it was shown in |11) the second boundary integral equation in (|17[) for the example problem (|49[l takes 
the form 

b(t)u(l,t) = b(t)v(l,t) + b(t) / —Wi{l,X,t-X)dX 

Jo M (52) 



= b(t)v(l,t) -b(t) K[t - X)g(X)dX + b(t) / K(t — X)b{X)u{l,X)dX, 
Jo Jo 

with 

oo 

K(t) = 2Y / e' [7r{2n ~ 1)/2]2t . (53) 

n=l 

The first equation in ([T7| is presented as follows 



ith 



u(x,t) =v(x,t) - / K r (t- X,x)g(X)dX+ / Ki(t - A, x)b(X)u(l, X)dX, (54) 



Kx(t, x) = 2j2(-l) n+2 e- [7r{2n - 1)/2]2t sin(-(2n - l)x), 



v(x, t) = 2 V e-M 2 "- 1 )/ 2 ! 2 ' sin( J(2n - l)x) f 
n=l 1/0 



uo(Osm(|(2n-l)Ode. 



(55) 



Remark 4.1 iVoie i/iai in i/ws particular case we can represent the integrand through the kernel K\(t — X,x) 
( obviously that K{t — X) = Ki(t — X, 1) ) analytically. In general case one can use the exponentially convergent 
algorithm for the operator exponential in (121|) like the ones from JSjj, fZlF- 
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Let us consider the particular case of the problem when 

7T 2 7T- 2 7T 

b(t)=c-—\ g(t)=e- 3 —\ u (ar) = sin(-x), 

with exact solution 

7T 2 7T 

t) = e~~* sin(— a;). 
Further we use to the equations ([521) . ([5"4"f the collocation method (12"2"1) . 

The results of computations are presented in the tables (HJ)-©. The first column indicates the collocation 
point tk, in the second column there are the errors of approximation in the boundary z x .k (i.e. for u(l,i)), 
in the third column there are the errors of approximation in the collocation points for x — ^. 



Point t 


£l 


£2 


.8535533905 


.28334234e-2 


.19383315e-2 


.1464466094 


.139098462e-l 


.75016794e-2 


Table 1: The error in the case n = 2, T = 1 


Point t 


El 


£2 


0.9619397662 


.36662211e-4 


.23005566e-4 


0.6913417161 


.34443339e-4 


.33521073e-4 


0.3086582838 


.42572982e-3 


.25374395e-3 


0.0380602337 


.23840042e-3 


.22946416e-4 


Table 2: The error in the case n = 4, T = 1 


Point t 




£2 


0.9903926402 


.46943895e-9 


.14327614e-9 


0.9157348061 


.14308953e-9 


.22763917e-9 


0.7777851165 


.1564038e-8 


.45939028e-9 


0.5975451610 


.2439358e-8 


.13481294e-9 


0.4024548389 


.9823893e-8 


.28128559e-8 


0.2222148834 


.1794445e-7 


.27310794e-8 


0.0842651938 


.3218373e-7 


.61284010e-8 


0.0096073597 


.1009601e-7 


.30125449e-ll 



Table 3: The error in the case n = 8, T = 1 
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